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Abstract 

In a general fractional factorial design, the n-levels of a factor are coded by the n-th 
roots of the unity. This device allows a full generalization to mixed-level designs of 
the theory of the polynomial indicator function which has already been introduced 
for two level designs by Fontana and the Authors (2000). the properties of orthogonal 
arrays and regular fractions are discussed. 
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1 Introduction 



Algebraic and geometric methods are widely used in the theory of the de- 
sign of experiments. A variety of these methods exist: real linear algebra, Z p 
ar ithme t ic, Ga lois Fields GF(p s ) arithmetic, where p is a prime number as 



in 



Dev and Mnkerjeel (119991 ) and lWu and Hamadal (1200(1 ). 



Bosd (|1947f). Se e , e.g., iRaktoe et al.l (|1981l) and the more recent books by 



Compl ex coding of l e vels has been used by many a u thors in various co n texts, 



f ll995h . lCollombierl (11999 ) and Ku and Wul (120011 ) 



see e.g . iBailevi ( 19821). iKobilinskv and Monodl (119911) , lEdmondsonl (119941 ) , iKobilinsky and Monod 
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The use of a new background, calle d Commutative Al g ebra or Polynomial 
Ring Algebra, was first advocated bvlPistone and Wynnl (119961 ) and later dis- 
cussed in detail i n iPistone et al.l (|2001l). O t her re leva nt general refer e nces are 
Robbianol (jl998l ). iRobbiano and Rogantinl (Il998l ) and lQaletto et all J20031). 



In the present paper, mixed- level (or asymmetric) designs with replicates are 
consid ered and the approach to th e two- level designs discussed in Font ana et al. 
f ll997h and iFontana et all (l2000h is generalized. In the latter, the fractional 
factorial design w as encoded in its indica tor function with respect to the full 
factorial design. In lTang and Deng! (11999I ). entities related to coefficients of the 
polynomial indicator function were independently introduced into the con- 
struction of a generalized wor d leng t h pat tern. The coefficients themselves 
were called J-characteristics in iTangl (1200 ll ). where it was shown that a two- 
level fractional design is uniquely determined by its J-characteristics. The rep- 
resentation of a fraction by its i ndicat or polynomial function was generalized 
to designs with replicates in |YeJ (120031 ) and extended to non t wo-level factors 
using orthogonal polynomials with an integer coding of levels in lCheng and Ye 
(l2004h . 



Sections 2 and 3 are a self-contained introduction of the indicator function 
representation of a factorial design using complex coding. The main results 
are in Sections 4 to 6. The properties of the indicator polynomial are dis- 
cussed in Section 4. If the factor levels are coded with the n-th roots of the 
unity, the coefficients of the indicator polynomial are related to many inter- 
esting properties of the fraction in a simple way: orthogonality among the 
factors and interactions, projectivity, aberration and regularity. Combinato- 
rial orthogonality vs. geometrical orthogonality is discussed in Section 5. A 
type of generalized regular fraction is defined and discussed in Section 6. The 
usual definition, where the number of levels is prime for all factors is extended 
to asymmetric design with any number of levels. With such a definition, all the 
monomial terms of any order are either orthogonal or totally aliased. However, 
our framework does not include the GF(jo s ) case. Some examples are shown 
in Section 7. 



A first partial draft of the present paper was presented in the GROSTAT 
V 2003 Workshop. S ome o f the results of Proposition 5 have been obtained 
independently by|YgJ (120041 ) . 



2 Coding of factor levels 



Let m be the number of factors of a design. We denote the factors by Aj, 
j = 1, ... ,m, and the number of levels of the factor Aj by rij. We consider 
only qualitative factors. 
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We denote the full factorial design by T>, V 
of all real responses defined on T> by TZ(T>). 



Ai X • • • x A m , and the space 



In some cases, it is of interest to code qualitative factors with numbers, es- 
pecially when the levels are ordered. Classical examples of numerical coding 



with rational numbers au £ Q are: (1) 



(2i 



rij 



l)/2 for odd rij and Qj 



i, or (2) a,ij 
2i — rij 



1, or (3) 



1 for even rij, see 

(jRaktoe et al.l . Il98ll . Tab. 4.1). The second case, where the coding takes value 
in the additive group Z n ., i.e. integers mod rij, is of special importance. We 
can define the important notion of regular fraction in such a coding. The third 
coding is the result of the orthogonalization of the linear term in the second 
coding with respect to the constant term. The coding — 1,+1 for two-level 
factors has a further property in that the yalues — 1,-1-1 form a mul tiplica - 
tive group. This prope rty w as widely us ed in iFontana et al.l (120001 ) , |Y§J (120031 ) , 
Tang and Pern] jl999h and|Tang] (l200lh . 



In the present paper, an approach is taked to parallel our theory for two-level 
factors with coding —1, +1. The n levels of a factor are coded by the complex 
solutions of the equation ( n = 1: 



2tt 

ujk = exp l l — k 

n 



k = 0, . . . , n — 1 . 



We denote such a factor with n levels by fl n , Q n = {u , . . . , cu n _i}. With such 
a coding, a complex orthonormal basis of the responses on the full factorial 
design is formed by all the monomials. 

For a basic reference to the algebra of the comple x field C an d of the n-th 
complex roots of the unity references can be made to lLan gj (119651 ); some useful 
points are collected in Section 8 below. 

As a = (3 mod n implies = , it is useful to introduce the residue class 
ring Z n and the notation [k] n for the residue of k mod n. For integer a, we 
obtain {oJk) a = V[ a k] n - The mapping 



Z n < — ► 



n n c c 



with 



k 



(2) 



is a group isomorphism of the additive group of Z n on the multiplicative group 
Q n C C. In other words, 

We drop the sub-n notation when there is no ambiguity. 
We denote by: 

• the number of points of the full factorial design, j^T> = IlfLi^j- 



L: the full factorial design with integer coding {0, 



,rij 



3 



and T> the full factorial design with complex coding: 

L = Z ni x • • • x Z nm and V = V 1 x • • • Vj ■ ■ ■ x V m with £>, = Q nj 

According to map (2), L is both the integer coded design and the exponent 
set of the complex coded design; 

• a, (3, . . . : the elements of L: 

L = {a = («i, . . . , a m ) : atj = 0, . . . , rij - 1, j = 1, . . . , m} ; 

that is, a is both a treatment combination in the integer coding and a 
multi-exponent of an interaction term; 

• [a - 0\: the m-tuple ([a x - A] ni , . . . , [atj - f3j] n , , . . . , [a m - f3 m ] nm ) ! the 
computation of the j-th element is in the ring Z n . . 



3 Responses on the design 



The responses on the design and the linear models are discussed in this section. 



Accor ding to the generalization of the algebraic approach by iFontana et al 



( 120001 ). the design T> is identified as the zero-set of the system of polynomial 
equations 

C"1 = , j = l,...,m. 
A complex response / on the design D is a C-valued function defined on T>. 
This response can be considered as the restriction to P of a complex polyno- 
mial. 

We denote by: 

• Xi] the z-th component function, which maps a point to its z-th component: 

X r . P 9 (&,..., Cm) I— ► Ci • 

The function Xi is called simple term or, by abuse of terminology, factor. 

• X a , with a G L: the interaction term Jf 1 ■ ■ • A"^ m , i.e. the function 

x a -. 2>9(k,...,Cm) ^ Cr---C m , <*eL. 

The function X a is a special response that we call monomial response or 
interaction term, in analogy with current terminology. 

In the following, we shall use the word term to indicate either a simple term 
or an interaction term. 

We say term X a has order (or order of interaction) k if k factors are involved, 
i.e. if the m-tuple a has k non-null entries. 
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If / is a response defined on V then its mean value on V, denoted by Ez>(f), 
is: 

E vU) = E /(O • 

We say that a response / is centered if Ex>(f) = 0. Two responses / and g are 
orthogonal on T> if Ex>(f ~g) = 0. 

it should be noticed that the set of all the responses is a complex Hilbert space 
with the Hermitian product f ■ g = Ex>(f g)- 

Two basic properties connect the algebra to the Hilbert structure, namely 

(1) x a W = x^-ft; 

(2) E V (X°) = 1, and E v (X a ) = for a ^ 0, see Section 8 Item (3). 

The set of functions {X a , a G L} is an orthonormal basis of the complex 
responses on design D. From properties (1) and (2) above it follows that: 



E v {X a XP) = E V {X 



1 if a = /3 
if a ^ (3 



Moreover, #L = #£>. 



Each response / can therefore be represented as a unique C-linear combination 
of constant, simple and interaction terms: 

f=J20 a X a , 9 a eC (3) 
where the coefficients are uniquely defined by: 6 a = Ed (^fX a ^j. In fact, 

E f (0x^(0 = E E 0px^(c) = E h E ^(0^(0 = e a . 

We can observe that a function is centered on T> if, and only if, 6 = 0. 
As # Q = E v {jX°\ , the conjugate of response / has the representation: 



/(C) = E * a (0 = E #[-^ a (0 • 

A response / is real valued if, and only if, 9 a = 9\- a ] f° r a h ct E L. 

We suggest the use of the roots of the unity because of the mathematical 
convenience we are going to show. In most of the applications, we are inter- 
ested in real valued responses, e.g. measurements, on the design points. Both 
the real vector space TZ(T>) and the complex vector space C(V) of the re- 
sponses on the design T> have a real basis, see (iKobilinskyl . Il990l . Prop. 3.1) 
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and iPistone and Rogantinl (120051 ) , where a special real basis that is common 
to both spaces is computed. The existence of a real basis implies the existence 
of real linear models even though the levels are complex. 



4 Fractions 



A fraction T is a subset of the design, T C T>. We can algebraically describe 
a fraction in two ways, namely using generating equations or the indicator 
polynomial function. 



4-1 Generating equations 



All fractions can be obtained by adding further polynomial equations, called 
generating equations, to the design equations Xj 3 — 1 = 0, for j — 1, . . . , m, 
in order to restrict the number of solutions. 



For e xample, let us consider a classical 3f n 2 regular fraction, see (jWu and Hamada 



20001 . Table 5A.1), coded with complex numbers according to the map in 
Equation (2). This fraction is defined by Xj —1 = for j = 1, . . . , 4, together 
with the generating equations X 1 X 2 X| = 1 and X 1 X|X 4 = 1. Such a rep- 
resentation of the fraction is classically termed "multiplicative" notation. In 
our approach, it is not a question of notation or formalism, but rather the 
equations are actually defined on the complex field C. As the recoding is a 
homomorphism from the additive group Z3 to the multiplicative group of C, 
then the additive generating equations in Z3 (of the form A + B + 2C = 
mod 3 and A + IB + D = mod 3) are mapped to the multiplicative equa- 
tions in C. In this case, the generating equations are binomial, i.e. polynomial 
with two terms. 

In the following, we consider general subsets of the full factorial design and, 
as a consequence, no special form of the generating equations is assumed. 



4-2 Responses defined on the fraction, indicator and counting functions 



The indicator polynomial was first introduc ed in Font ana et al.l (119971 ) to de- 
scribe a fraction. In the two-level case, |YgJ (120031 ) suggested generalizing the 
idea of indicator function to fractions with replicates. However, the single repli- 
cate case has special features, mainly because, in such the equivalent 
description with generating equations is available. For coherence with gen- 
eral mathematical terminology, we have maintained the indicator name, and 
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introduced the new name, that is, counting function for the replicate case. 
The design with replicates associated to a counting function can be consid- 
ered a multi-subset T of the design T>, or an array with repeated rows. In the 
following, we also use the name "fraction" in this extended sense. 

Definition 1 (Indicator function and counting function) The counting 
function R of a fraction T is a response defined on T> so that for each (eP, 
R(C) equals the number of appearances of £ in the fraction. 

A 0-1 valued counting function is called indicator function F of a single repli- 
cate fraction T . 

We denote the coefficients of the representation of R onT> using the monomial 
basis by b a : 

r(C) = J2 b » xa (0 C € v . 



A polynomial function R is a counting function of some fraction T with repli- 
cates up to r if, and only if, R[R — 1) ■ • • (R — r) = on D. In particular a 
function F is an indicator function if, and only if, F 2 — F = on T>. 

If F is the indicator function of the fraction JF, F — 1 = is a set of generating 
equations of the same fraction. 

As the counting function is real valued, we obtain b a = 6[_ a j. 

If / is a response on T> then its mean value on !F, denoted by Ep(f), is: 

EAf) = ^ E /(O = || MR f) 

where j^T is the total number of treatment combinations of the fraction, 

Proposition 1 (1) The coefficients b a of the the counting function of a frac- 
tion T are: 

in particular, b is the ratio between the number of points of the fraction 
and those of the design. 
(2) In a single replicate fraction, the coefficients b a of the indicator function 
are related according to: 

b a = ^2 bp b[ a _p] . 
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(3) If T and J 7 ' are complementary fractions without replications andb a and 
b' a are the coefficients of the respective indicator functions, bp, = 1 — b' Q 
and b a = —b' a . 

PROOF. Item (1) follows from : 

E ^CJ E R ^CJ = E E b^{Qx^) J2b a = #vb a . 

Item (2) follows from relation F = F 2 . In fact: 

e b a x« = y: bpx? y: b^ = y: b^x^ = 

= E E 6M ft = EE^Mi^. 

a \P+l\=<x a & 

Item (3) follows from F' = 1 — F. □ 



4-3 Orthogonal responses on a fraction 



In this section, we discuss the general case of fractions JF with or without 
replicates. As in the full design case, we say that a response / is centered on 
a fraction T if Ej?(f) = Ex>{R f) = and we say that two responses / and 
g are orthogonal on T if Ej?(f ~g) = E-p(R f g) = 0, i.e. the response / ~g is 
centered. 

It should be noticed that the term "orthogonal" refers to vector orthogonality 
with respect to a given Hermitian product. The standard practise in orthogo- 
nal array literature, however, is to define an array as orthogonal when all the 
l evel combinations ap pear equally often in relevant subsets of columns, e.g. 



( iHedayat et al.l . ll999l . Def. 1.1). Vector orthogonality is affected by the coding 
of the levels, while the definition of orthogonal array is purely combinatorial. A 
characterization of orthogonal arrays can be based on vector orthogonality of 
special responses. This section and the next one are devoted to discussing how 
the choice of complex coding makes such a characterization as straightforward 
as in the classical two-level case with coding -1,+1 . 

Proposition 2 Let R = J2aeLb a X a be the counting function of a fraction T . 

(1) The term X a is centered on T if, and only if, b a = &r_ a i = 0. 

(2) The terms X a and X^ are orthogonal on T if, and only if, b\ a -p\ = 0; 

(3) If X a is centered then, for each (3 and 7 such that a = [(3 — 7] or a = 
[7-/5], X 13 is orthogonal to X 1 . 

(4) A fraction T is self- conjugate, that is, R(() = R{() for any ( £ T>, if, 
and only if, the coefficients b a are real for all a £ L. 



S 



Proof. The first three Items follow easily from Proposition 1. 



For the Item (4), we obtain: 

R(0 = E b a x°(C) = E h- a ]X ha] (0 = E » [ ~ a] (C) 
R(0 = E b a X a (() = E b a x [ ~ a \c) . 

Therefore R(() = -R(C) if; an d only if, b a = b a . It should be noticed that the 
same applies to all real valued responses. □ 

Interest in self-conjugate fractions concerns th e existence o f a rea l valued linear 
basis of the response space, as explained in dKobilinskvl . Il99(l Prop. 3.1). It 



follows that it is possible to fit a real linear model on such a fraction, even 
though the levels have complex coding. 

An important property of the centered responses follows from the structure 
of the roots of the unity as a cyclical group. This connects the combinatorial 
properties to the coefficients b a 7 s through the following two basic properties 
which hold true for the full design T>. 

P-l Let Xi be a simple term with level set Q n . Let us define s = n/gcd(r,n) 

and let Q s be the set of the s-th roots of the unity. The term XJ takes 

all the values of fl s equally often. 
P-2 Let X a = X: 31 • ■ • X 4 Jk be an interaction term of order k where X, H 

takes values in Vt s . Let us define s = \cm{sj 1 , . . . , Sj k }. The term X a 

takes values in fl s equally often. 

Let X a be a term with level set Q s on the design T>. Let rj, be the number of 
times X a takes the value u>k on J 7 , k — 0, . . . , s — 1. The polynomial P(C) is 
associated to the sequence (rk)k=o,...,s-i so that: 

nO = X>C fc withCeC. 

fc=0 

It should be noticed that 



See iLangl (119651 ) and the Appendix for a review of the properties of such a 



polynomial P. 

Proposition 3 Let X a be a term with level set Q s on full design T>. 
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(1) X a is centered on T if, and only if, 

P(C) = *.(C)*(0 

where $ s is the cyclotomic polynomial of the s-roots of the unity and 
is a suitable polynomial with integer coefficients. 

(2) Let s be prime. Therefore, the term X a is centered on T if, and only if, 
its s levels appear equally often: 

r = ■ ■ • = r s _i = r 

(3) Let s = p^ 1 pj d , with p^ prime, for i = l,...,d. The term X a 

is centered on T if, and only if, the following equivalent conditions are 
satisfied. 

(a) The remainder 

H({) = P(C) mod $ S (C) , 

whose coefficients are integer combination of r^, k = 0, . . . ,s — 1, is 
identically zero. 

(b ) The polynomial of degree s 

HO = P«)U*d(C) mod(C s -l), 

d\s 

whose coefficients are integer combination of the replicates r^, k = 
0, . . . , s — 1, is identically zero. The indices of the product are the d's 
that divide s. 

(4) Let gi be an indicator of a subgroup or of a lateral of a subgroup of Q s ; 
i.e.: g { = (g a , . . . ,gij, . . . ,g is ), g^ G {0, 1} ; such that {k : g ik = 1} is a 
subgroup or a lateral of a subgroup of Q s . 

If the vector of level replicates (ro,7*i, . . . , r s _i) is a combination with 
positive weights of g^: 

(r ,ri, . . . ,r a _i) = ^ a * 9i withoi^N 

X a is centered. 

Proof. (1) As Uk = oj\, the assumption J2k r k^k = is equivalent to 
P(o>i) = 0. From Section 8, Items 4 and 5, we know that this implies 
that P(uj) = for all primitive s-roots of the unity, that is, P(C) is 
divisible by the cyclotomic polynomial $ s . 

(2) If s is a prime number, the cyclotomic polynomial is $ S (C) = J2k^oC k - 
The polynomial P(C) is divided by the cyclotomic polynomial, and P{() 
and $ S (C) have the same degree, therefore r s _i > and P(C) = r s-i$((), 
so that 7*o = • • • = r s _i. 

(3) The divisibility shown in Item 1 is equivalent to the condition of null 
remainder. Such a remainder is easily computed as the reduction of the 
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polynomial P(C) mod $ S (C)- According to the same condition and Equa- 
tion (8), we obtain that -P(C) is divisible by ( s — 1, therefore it also equals 
mod ( s — 1. 

(4) If Q p is a prime subgroup of Q s , then J2wen p ^ = 0. Now let us assume that 
the replicates on a primitive subgroup Q Pi are 1. Therefore J2u&n . u — 
according the equation in Item (3). The same occurs in the case of the 
laterals and the sum of such cases. 

□ 

Example 

Let us consider the case s = 6. This situation occurs in the case of mixed-level 
factorial designs with both three-level factors and two-level factors. In this 
case, the cyclotomic polynomial is $6(0 = C 2 — C + 1 whose roots are uj\ and 
005. The remainder is 

#(0 = I>C fc mod$ 6 (C) 

fc=0 

= r + riC + r 2 ( 2 + r 3 ( 3 + r 4 ( 4 + r 5 ( 5 mod (( 2 - C + l) 
= (n + r 2 -r 4 - r 5 )C + (r - r 2 - r 3 + r 5 ) 

The condition H(Q = implies the following relations concerning the numbers 
of replicates: r +ri = r 3 +r 4 , ri+r 2 = r i +r 5 , r 2 + r 3 = r + r 5 , where the 
first one follows by summing of the second with the third one. Equivalently: 

r - r 3 = r 4 - n = r 2 - r 5 . (4) 

Let us consider the replicates corresponding to the sub-group {a>o, ^2, 0J4} and 
denote the min{ro, r 2 , r 4 } by mi. We then consider the replicates correspond- 
ing to the lateral of the previous sub-group {uj\, L03, u$} and we denote by m 2 
the min{r 1; r 3 , r 5 }. We consider the new vector of the replicates: 

= (r - mi, r\ - m 2 , r 2 - m 1 , r 3 - m 2 , r 4 - mi, r 5 - m 2 ) 
= r - 774(1, 0, 1, 0, 1, 0) - m 2 (0, 1, 0, 1, 0, 1) 

The vector r' satisfies Equation (4). 

As at least Tq, r' 2 or r 4 is zero, the common value in Equation (4) is zero 
or negative. Moreover, as at least r[, r 3 or r' 5 is zero, the common value in 
Equations (4) is zero or positive. The common value is therefore zero and 
r = r 3' r 'l = r 4' r> 2 = r 5 an d 

r' = r' (l, 0, 0, 1, 0, 0) + r{(0, 1, 0, 0, 1, 0) + r 2 (0, 0, 1, 0, 0, 1) 
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A term is therefore centered if the vector of the replicates is of the form: 

(r , . . . , r B ) = ai(l, 0, 0, 1, 0, 0) + a 2 (0, 1, 0, 0, 1, 0) 

+ a 3 (0, 0, 1, 0, 0, 1) + a 4 (l, 0, 1, 0, 1, 0) + a 5 (0, 1, 0, 1, 0, 1) 

with dj non negative integers. There are 5 generating integer vectors of the 
replicate vector. 

It should be noticed that if the number of levels of X a is not prime, Ejr(X a ) = 
does not imply Ejr[X ra ) = 0. In the previous six-level example, if X a is 
centered, the vector of replicates of X 2a is of the form (2ai + a 4 + a^, 0, 2a 2 + 
a 4 + a 5 , 0, 2a 3 + a 4 + a 5 , 0) and X 2a is centered only if a\ = a 2 = a 3 . 



5 Orthogonal arrays 



In this sectionn we discuss the relations between the coefficients b a , a £ L, of 
the counting function and the property of being an orthogonal array. Let 

OA(n, S f,..., S £T,t) 

be a mixed-level orthogonal array with n rows and m columns, m = p\ + • ■ • + 
p m , in which pi columns have s\ sy mbols Pk columns h ave s m symbols, and 



with strength t, as defined e.g. in (IWu and Hamadal . |2000| . p. 260). Strength t 



means that, for any t columns of the matrix design, all possible combinations 
of symbols appear equally often in the matrix. 

Definition 2 Let I be a non-empty subset of {l,...,m}, and let J be its 

complement set, J = I c . Let T>i and T>j be the corresponding full factorial 
designs over the I -factors and the J -factors, so that T> = T>i xT>j. Let T be 
a fraction of T> and let Ti and Tj be its projections. 

(1) A fraction T factorially projects on the I -factors if Ti = s Vj, that is, 
the projection is a full factorial design where each point appears s times. 

(2) A fraction T is a mixed orthogonal array of strength t if it factorially 
projects on any I -factors with #1 = t. 

Using the notations of Definition 2, for each point £ of a complex coded fraction 
JF, we consider the decomposition £ = (j) and we denote the counting 
function restricted to the /-factors of a fraction by Rj, i.e. Ri((i) is the number 
of points in T whose projection on the /-factors is (j. 

We denote the sub-set of the exponents restricted to the /-factors by Lj and 
an element of Lj by otf. 

Li = {ai = (ai, . . . , aij, . . . , a m ), a,- = if j e J} . 
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Therefore, for each a E L and ( 6 P: a = ar+aj and X a (() = X ai (Ci)X aj (Cj). 
We denote the cardinalities of the projected designs by and #T>j. 

Proposition 4 

(1) The number of replicates of the points of a fraction projected onto the 
I-factors is: 

(2) A fraction factorially projects onto the /-factors if, and only if, 

Ri(Ci) = #Vj b = ^- for all (i • 

This is equivalent to all the coefficients of the counting function involving 
only the I-factors being 0: 

b ai = with ai G Li, at/ ^ (0, 0, . . . , 0) . 

In such a case, the levels of a factor Xi, i G / , appear equally often in T . 

(3) If there exists a subset J of {1, . . . , m} such that the J -factors appear in 
all the non null elements of the counting function, the fraction factorially 
projects onto the /-factors, with I = J c . 

(4) A fraction is an orthogonal array of strength t if, and only if, all the 
coefficients of the counting function up to the order t are zero: 

b a = V a of order up to t, a ^ (0,0,...,0) . 

Proof. (1) We obtain: 

/MC/)= E tf(C/,G)= E E^ Q (C/,G) 

= E E^ Q/ (C/P^(G) 

= E b ai x^(o) e * QJ (C)+ E b a x a ^o e x aj (0 ■ 

The thesis follows from E(j£Vj X aj {(j) = if aj ^ (0, 0, . . . , 0) and 
E cje vj X aj (G) = if olj = (0, 0, . . . , 0). 

(2) The number of replicates of the points of the fraction projected onto the 
/-factors, Ri((i) = ^PjEa/^a/ X a '(d), is a polynomial and it is a 
constant if all the coefficients b an with a>i ^ (0, 0, . . . , 0), are zero. 

(3) This condition implies that the b a /s are zero, if ai ^ (0,0, .. . ,0), and 
the thesis follows from the previous item. 

(4) This item follows from the previous items and the definition. 

□ 
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Remarks 

(1) If a fraction factorially projects onto the /-factors, its cardinality must 
be equal to, or a multiple of the cardinality of Vj. 

(2) If the number of levels of each factors is a prime, the condition b a . = 
for each i 6 J and < ati < rii — 1 in Items (2) and (3) of the previous 
Proposition, simplify to [Xj) = 0, according to Item (2) of Proposition 
3. 



6 Regular fractions: a partial generalization to mixed-level design 



A short review of the theory of regular fractions is here made from the view 
point of the present paper. Vario us definitions of reg ular fractio n appear in 



litera ture, e.g. in the boo k s by ( Raktoe et al. . 1981, p. 123 ) , (Collombieri 



f Wu and Hamadal . 12000 



19961. p. 125). flKobilinskvl . Il997t p. 70). flDev and Mukerieel . Il999l . p. 164) 



p. 305). To our knowledge, all the definitions are 
known to be equivalent if all the factors have the same prime number of lev- 
els, n = p. The definition based on Galois Field computations is given for 
n = p s power of a prime number. All definitions assume symmetric factorial 
design, i.e. all the factors have the same number of levels. 



Regular fraction designs are usually considered for qualitative factors, where 
the coding of the levels is arbitrary. The integer cod ing, th e GF(p s ) coding, 
and the roots of the unity coding, as introduced by iBaileyl (119821 ) and used 
extensively in this paper, can all be used. Each of those codings is associated 
to specific ways of characterizing a fraction, and even more important for 
us, to a specific basis for the responses. One of the possible definitions of a 
regular fraction refers to the property of non-existence of partial confounding 
of simple and interaction terms, and this property has to be asso ciated to a 
specific basis, as explicitly pointed out in IWu and Hamadal (120001 ) . 



In our approach, we use polynomial algebra with complex coefficients, the n- 
roots of the unit coding, and the idea of indicator polynomial function, and 
we make no assumption about the number of levels. In the specific coding we 
use, the indicator polynomial is actually a linear combination of monomial 
terms which are centered and orthogonal on the full factorial design. We refer 
to such a basis to state the no-partial confounding property. 

The definition of the regular fraction is hereafter generalized in the symmetric 
case with a prime number of levels. The new setting includes asymmetric 
design with any number of levels. Proposition 5 below does not include regular 
fractions defined in GF(p s ). A full discussion of this point shall be published 
elsewhere. 
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We consider a fraction without replicates. Let n = lcmjri!, . . . ,n m }. It should 
be recalled that Q n is the set of the n-th roots of the unity, Q n = {uq, . . . , u n _i}. 
Let £ be a subset of exponents, £ C L = 1\ x ■ • • x Z m , containing (0, . . . , 0) 
and let I be its cardinality (/ > 0). Let e be a map from £ to fl n , e : £ — > fl n . 

Definition 3 A fraction T is regular if 

(1) £ is a sub-group of L, 

(2) e is a group homomorphism, e([a + [3]) = e(a) e(/3) for each 

(3) the equations 

X a = e(a) a G £ (5) 

define the fraction T , i.e. they are a set of generating equations, according 
to Section J^.l. Equations (5) are also called the defining equations of T . 

IfTi is a minimal generator of the group £, Equations X a = e(a), a G Ti C £, 
are called a minimal set of generating equations. 

It should be noticed that our situation is general because the values e(a) can 
be different from 1. From items (1) and (2) it follows that a necessary condition 
is that the e(a)'s must belong to the subgroup spanned by the values of X a . 
For example, for n\ = n2 = n = 6, an equation such as XfX% = 002 cannot be 
a defining equation. 

For example, in the fraction of Section 4.1, we have: H. — {(1, 1, 2, 0), (1, 2, 0, 1)} 
and e(l, 1, 2, 0) = e(l, 2, 0, 1) = uj q = 1. The set £ is: {(0, 0, 0, 0), (0, 1, 1, 2), 
(0, 2, 2, 1), (1, 1, 2, 0), (2, 2, 1, 0), (1, 2, 0, 1), (2, 1, 0, 2), (1, 0, 1, 1), (2, 0, 2, 2)}. 

Proposition 5 Let T be a fraction. The following statements are equivalent: 

(1) Fraction T is regular according to Definition 3. 

(2) The indicator function of the fraction has the form 

^(C) = )E^)^(0 

where £ is a given subset of L and e : £ — » Q n is a given mapping. 

(3) For each G L, the parametric functions represented on JF by the 
terms X a and X 13 are either orthogonal or totally confounded. 

Proof. First we prove the equivalence between (1) and (2). 
(1) =► (2). 

Let T be a regular fraction and let X a = e(a) be its defining equations with 
a G £, £ a sub-group of L and e a homomorphism. 
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If, and only if, ( E T: 



= E l* a (0 - e(«)| 2 = £ (X°(C) - e(a)) (X«(0 - e(a)) 



£ (X a (OX«(() + e(a)e(a) - e(a)X<*{() - e(a)X a (() 



2 / - E e(a)X-(C) - E *°(C) = 2 f / - E ^(o) 



therefore 

1 



E e («) X °(0 -1 = if, and only if, ( E T . 



The function F = y X) ag £ e(a) I" is an indicator function, as it can be shown 
that F = F 2 on V. In fact, C is a sub-group of L and e is a homomorphism; 
therefore: 

^ 2 = 4e E^M* [a+/3] = 4 EE*+^])^ = 



- E « <l) x^ = f 

1 ~f€C 



It follows that F is the indicator function of J 7 , and b a = j^-, for all a E C 
(2) =► (1). 

It should be noticed that an indicator function is real valued, therefore F = F. 

7El^(0-e(a)| 2 = 2-F(C)-F(C) = 2-2F(0 = (° ° U £ 

Equations X a = e(a), with a G £, define the fraction JF as the generating 
equations of a regular fraction. It is easy to see that £ is a group. In fact, if 7 = 
[a + 0\£ C, there exists one ( such that X~<(() = X a (()X^(C) = e(a)e(/3) C 
Q n and the value e(a)e(/3) only depends on 7. By repeating the previous proof, 
the uniqueness of the polynomial representation of the indicator function leads 
a contradiction. 

Now we prove the equivalence between (2) and (3). 

(2) =► (3) 

The non-zero coefficients of the indicator function are of the form b a = e(a)/l. 

We consider two terms X a and X 13 with a, j3 G L. If [a — /3] E* C then X a and 
X? are orthogonal on T as the coefficient b[ a -p] of the indicator function equals 
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0. If [a — (3] G £ then X a and X 13 are confounded because X^- a ® = e([a — j3]); 
therefore X a = e([a - 0\) X 13 . 

(3) =>■ (2). 

Let £ be the set of exponents of the terms confounded with a constant: 
£ = {a G L : X a = constant = e(a), e(a) G f2 n } . 



For each a G £, 6 a = e(o;) 6 - For each a ^ £, because of the assumption, X a 
is orthogonal to X°, therefore b a = 0. □ 

Corollary 1 Lei J 7 be a regular fraction with X a = 1 for all the defining 
equations. T is therefore self- conjugate and a multiplicative subgroup ofV. 

Proof. It follows from Prop. 2 Item 4. □ 



The following proposition extends a result presented in iFontana et al. 
for the two level case. 

Proposition 6 Let T be a fraction with indicator function F. We denote the 
set of the exponents a such that ^ = e(a) G f2 n by £ . The indicator function 
can be written as 



F(0 = &oE e («) x °(0 + ^(0 C e , £ n K = . 

a££ /3e/C 

Ji follows that £ is a subgroup and the equations X a = e(a), wit/i a G £, 
are the defining equations of the smallest regular fraction T r containing T 
restricted to the factors involved in the £- exponents. 

Proof. The coefficients b a , a G £, of the indicator function F are of the form 
boe(oe). Therefore, from the extremality of n-th roots of the unity, X a (() = 
e(a) if C G T and X a (()F(Q = e(a)F(() for each ( G V and £ is a group. 

We denote the indicator function of T r by F r . For each (fDwe have: 

nom) = \no e * q (o = y e * q (o no = 

= yE^)e(«) F(C) = jjF(C). 

The relation F^F^) = F(£) implies F C F r . The fraction F r is minimal 
because we have collected all the terms confounded with a constant. □ 
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Remark 



Given generating equations X ai = 1, . . . , X ah = 1, with {a±, . . . , ah} = Ti. C 
Z m x • • • x Z„ m , the corresponding fraction JF is a subgroup of fl ni x • • • x fi nm . 
If the same fraction is represented in the additive notation, such a set of 
treatment combin ations is the principal b l ock of a single re p licate gene ralized 
cyclic design, see IJohn and Dean! (119751 ). iDean and John! (119751 ) and [Lewis 
( 119791 ). A complex vector of the form 



i 2tt^ 

e ™i , e 



l 2-rr — 

~>2 



i 2n^- 



with < k{ < n i} % — 1, 



m 



is in fact a solution of the generating equations if, and only if, 



Y^j=\ a ij Ij kj = mod s 



< kj < rij 



(6) 



with s = lcm{ni, . . . , n m } and jj = 



A set of generators can be computed from Equation (6). It should be noticed 
that the following equivalent i nteger linear pro gramming problem does not 
involve computation mod s, see ISchrijverl ( 120021 ) 



E^Li Oiij 7j kj - sq = 
< kj <nj , q > 



(7) 



In iLewisi ( 119821 ). the mon omial part of our de fining equations is called defin- 
ing contrast, according to iBailey et al.l (119771 ). The paper contains extensive 
tables of the generator subgroups of the treatment combinations and the cor- 
responding defining contrasts. 



Viceversa, given a set of generators of the treatment combinations, 

{&!, ...,b r \bi = i:L , . . -Mm)} , 

Equation (6) with indeterminates an 

YlJLi Ij hj oij = mod s 
< OLj < rij 



produces generating equations for the fraction. 
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7 Examples 



A regular fractions with n = 3. 

Let us consider the classical 3 4-2 fraction of Section 4.1. Its indicator function 
is: 



F — — (l + X 2 X 3 X^ + X\X^X\ + X±X 2 X^ + X\X\X 3 

-\-X\X\X± + XfX 2 X% + X1X3X4 + X^X^X^j . 

We can observe that the coefficients are all equal to g. The minimum order of 
interactions that appear in the indicator function is 3, therefore the fraction 
is an orthogonal array of strength 2. All the defining equations are of the form 
X a = 1, therefore the fraction is self-conjugate. 



A regular fraction with n = 6. 

Let us consider a 6 3 design. From property [P-2] of Section 4.3, the terms 
X a take values either in fl e or in one of the two subgroups either {l,dj 3 } or 
{1,^2,^4}- 

Let T be a fraction whose generating equations are: X\X\X\ = u 3 and 
X 4 X 2 4 Xf = u 2 . In this case we have: H = {(3, 3, 3), (4, 4, 2)} and e(3, 3, 3) = 
cu 3 , e(4, 4, 2) = u 2 . The set C is: {(0, 0, 0), (3, 3, 3), (4, 4, 2), (2, 2, 4), (1, 1, 5), (5, 1, 1)}. 
The full factorial design has 216 points and the fraction has 36 points. The 
indicator function is: 

F = - (l + ou 3 XfXlXi + u±X\X\Xl + uj 2 X\X\Xl + uj x X x X 2 X\ + ou 5 XfX 2 X 3 ) 
It should be noticed that this fraction is an OA(36, 6 3 , 2). 



An OA(18,2 1 3 7 ,2). 

We c onsider the fraction of a 2x3 7 design with 18 runs, taken from (IWu and Hamada . 



20001 . Table 7C.2) and recoded with complex levels. Here X\ takes values in 
Q 2 , Xi, with i = 2, ... ,8, and their interactions take values in f2 3 , and the 
interactions involving X± take values in Q e . 

All the 4374 X a terms of the fraction have been computed in SAS using 1, 2 , 
Z3 and 7jq arithmetic. The replicates of the values in the relevant have then 
been computed for each terms. We found: 

(1) 3303 centered responses. These are characterized by Proposition 3. The 
replicates are of the type: (9, 9), (6, 6, 6), (3, 3, 3, 3, 3, 3) and (9, 0, 0, 9, 0, 0). 
We have: 

(a) the two- level simple term and 1728 terms involving only the three- 
level factors (14 of order 1, 84 of order 2, 198 of order 3, 422 of order 
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4, 564 of order 5, 342 of order 6 and 104 of order 7); 
(b) 1574 terms involving both the two-level factor and the three-level 
factors (14 of order 2, 66 of order 3, 188 of order 4, 398 of order 5, 
492 of order 6, 324 of order 7 and 92 of order 8). 

(2) 9 terms with corresponding b a coefficients equal to b = ^f? = 3~ 5 ; 

(3) 1062 terms with corresponding coefficients different from zero and bo: 450 
terms involving only the three-level factors (80 of order 3, 138 of order 4, 
108 of order 5, 100 of order 6 and 24 of order 7) and 612 terms involving 
both the two-level factor and the three-level factors (18 of order 3, 92 of 
order 4, 162 of order 5, 180 of order 6, 124 of order 7 and 36 of order 8). 

Some statistical properties of the fraction are: 

(1) Analyzing the centered responses we can observe that: 

(a) All the 15 simple terms are centered. 

All the 98 interactions of order 2 (84 involving only the three-level 
factors and 14 also involving the two-level factor) are centered. This 
implies that both the "linear" terms and the "quadratic" terms of the 
three-level factors are mutually orthogonal and they are orthogonal 
to the two-level factor. 

The fraction is a mixed orthogonal array of strength 2. 

(b) The fraction factorially projects onto the following factor subsets: 



All the terms of order 1, 2 and 3 involving the same set of factors are 
in fact centered. 

(c) The minimal regular fraction containing our fraction restricted to the 
three-level factors has the following defining relations: 



(d) The non centered terms have levels in f2 6 and in Q 3 . 
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8 Appendix: Algebra of the n-th roots of the unity. 

We hereafter list some facts concerning the algebra of the complex n-th roots 
of the unity, for ease of reference. 

(1) The conjugate of a n-th root of the unity equals its inverse: uJk = ui^ 1 = 
uj\-k] for all cUfc G fi n . 

(2) If C + u m , we obtain: YlVo k*J( ~ "k) = £^ = EU^C*" 1 ^ 
where the last equality follows from algebraic computation. Therefore, 
for C = uj m : 

n-l 

and especially: ]lfe=i(l — <^fe) — n - 

(3) We have: C n -1 = (C~^o) ■ ■ • (C-"»-i) = T^{-l) n ~ k S n - k K • • • , w n -i) C h 
where SVi-fc (x , . . . , is the elementary symmetric polynomial of or- 
der n — k. We therefore obtain the following notable cases: 

- Si [L)o, . . . , CJ n _i) = J2k = 

- 1S2 (tUo, . . . , CJ n _i) = YliKm U m = 

- S n (u) , . . . , cj n _i) = Uk ^k = (-l) n+1 

where the indices of the sums and products are from to n — 1. 

(4) Let u be a primitive n-th root of the unity, that is, a generator of Q n as 
a cyclic group: {l,cu, u 2 , . . . , u n ~ 1 } = Q n . 

The root uj p G Q n is primitive if p is relatively prime with n. In particular, 
ui is a primitive root and, for Uk G f2 n , we obtain: Uk = (<^i) fc - If u is 
a prime number, all the roots of the unity, except 1, are primitive roots. 
The number of the primitive n-th roots of the unity is denoted by 0(n). 

(5) Given an algebraic number x, the unique irreducible monic polynomial 
of the smallest degree with rational coefficients P such that P(x) = 
and whose leading coefficient is 1, is called the minimal polynomial of x. 
The minimal polynomial of a primitive n-th root of the unity is called 
the cyclotomic polynomial $ n (C) an d its degree is <f>(ri): 

$n(C) = II(C — Up) , ( G C, uj p G fl n primitive n-th root of the unity. 

p 

If n is prime, the minimal polynomial of a primitive n-th root of the unity 



n-l 
h=0 
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is $„(£) = C™- 1 + ( n ~ 2 + • • • + 1. Moreover: 

C n -l = $ n (C)----$d(C)----$i(C) where d divides n . (8) 

(6) The recoding in Equation (2) is a polynomial function of degree n — 1 
and complex coefficients in both directions: 

E 1 Uh=o,h^s {x — h) 
uj s n _{ - ( tt, x = k G {0, . . . ,n - 1} 



s=0 



n— 1 n— 1 



n^o,^ s (« - *o 



= - £ h Y. s u [8-sh] , ( = uj k ett n . (9) 



71 h=0 8=1 



The last Equation follows from 

and from the properties of the n-th roots of the unity, see Item 2. 
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